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ABSTRACT 


Hie  first  temperature  derivatives  of  the  fundamental  elastic  constants 
of  quartz  are  employed  along  with  the  thermally-induced  biasing  strains  in 
the  equation  for  the  first  perturbation  of  the  eigenvalue  for  the  linear 
electroelastic  equations  for  small  fields  superposed  on  a  bias  to  calculate 
the  resulting  change  in  surface  wave  velocity  with  temperature.  Since  the 
description  employed  is  referred  to  a  fixed  reference  state,  the  geometry 
does  not  change  and  the  temperature  coefficient  of  velocity  is  the  negative 
of  the  temperature  coefficient  of  delay.  In  all  earlier  work  a  variable 
temperature  dependent  geometric  state  was  considered,  but  the  attendant 
shearing  or  skewing  of  the  coordinate  axes  was  not  properly  included.  The 
temperature  dependence  of  the  velocity  of  surface  waves  propagating  in 
various  directions  on  various  cuts  of  quartz  is  obtained.  For  cases  in 
which  measurements  are  available  the  calculations  are  shown  to  be  in  sub¬ 
stantially  better  agreement  with  the  measurements  than  previous  calculations. 


.Accession  For 

Niio  Grin 
Z  ‘.  i- 


Work  performed  while  at  Rensselaer  Polytechnic  Institute 


1.  Introduction 


1  2 

Previous  work  ’  on  the  temperature  dependence  of  the  velocity  of 
surface  waves  in  quartz  has  employed  the  temperature  derivatives  of  the 
effective  elastic  constants  of  quartz,  which  are  referred  to  the  variable 
temperature  dependent  intermediate  position  rather  than  the  fixed  reference 
position,  to  which  the  fundamental  elastic  constants  are  referred.  Hot  only 
does  the  intermediate  coordinate  system  extend  or  contract  under  a  tempera¬ 
ture  change,  but  it  shears  or  skews  as  well.  This  shearing  or  skewing  of 

3 

the  axes  was  omitted  in  the  original  determination  of  the  temperature 

derivatives  of  the  effective  elastic  constants.  Although,  in  the  existing 
1  2 

treatments  ’  of  the  temperature  dependence  of  surface  wave  velocity,  the 

change  in  density  and  path  length  were  properly  included,  the  shearing  or 

skewing  of  the  coordinate  axes  was  omitted. 

Recently,  the  first  temperature  derivatives  of  the  fundamental  elastic 

4  3 

constants  of  quartz  were  determined  from  the  original  data  from  which  the 

3 

temperature  derivatives  of  the  effective  constants  had  been  obtained  .  Also, 

a  perturbation  analysis  of  the  linear  electroelastic  equations  for  small 

5  6 

fields  superposed  on  a  bias  has  been  performed  and  the  equation  for  the 
first  perturbation  of  the  eigenvalue  has  been  obtained.  Hie  aforementioned 
first  temperature  derivatives  of  the  fundamental  elastic  constants  of  quartz 
are  employed  along  with  the  thermally-induced  biasing  strains  in  the  equation 
for  the  first  perturbation  of  the  eigenvalue^  to  calculate  the  resulting 
change  in  surface  wave  velocity  with  temperature. 

Since  the  description  employed  is  referred  to  the  fixed  reference  state 
at  the  reference  temperature,  the  mass  density  and  geometry  do  not  change  and 
the  temperature  coefficient  of  natural  velocity,  which  we  determine,  is  the 
negative  of  the  temperature  coefficient  of  delay.  Hie  temperature  dependence 


2 


of  the  actual,  in  addition  to  that  of  the  natural,  velocity  of  surface  waves 
on  a  number  of  cuts  of  quartz  is  obtained  and  compared  with  the  previous 
calculations  as  well  as  measured  values,  in  addition,  the  temperature 
coefficient  of  delay  and  the  associated  power  flow  angle  are  calculated  as 
a  function  of  propagation  direction  for  two  interesting  orientations  of 
quartz  substrates.  All  results  obtained  here  are  limited  to  linear  behavior 
with  temperature  because  only  the  first  temperature  derivatives  of  the 

4 

fundamental  elastic  constants  of  quartz  are  known  .  As  noted  in  an  earlier 
7 

work  the  surface  wave  velocity  can  also  be  obtained  simply  by  solving  the 
appropriate  linear  surface  wave  boundary  value  problem  in  the  usual  manner 
at  both  the  reference  temperature  and  present  temperature  when  the  tempera¬ 
ture  dependent  constants  have  been  determined  at  the  present  temperature. 

The  perturbation  procedure  is  being  employed  here  because  it  enables  the 
change  in  surface  wave  velocity  to  be  calculated  directly,  thereby  resulting 
in  greater  accuracy  for  the  same  number  of  significant  figures  and  a  sub¬ 
stantial  reduction  in  computer  time. 

2.  Linear  Electroelastic  Equations  for  Small  Fields  Superposed 
on  a  Thermally  Induced  Bias 

Die  linear  electroelastic  equations  for  small  fields  superposed  on  a 
bias  may  be  written  in  the  form^’® 

V,i'p°V 

where,  since  only  mechanical  nonlinearities  are  of  interest  here, 

*I/y  ”  5l^MVUV,  M  +  eMLy^,  M  * 

V  eiMVUV,M~  eU^,  MJ 


(2.2) 


3. 


lIiyMV  "  §iiyMV  +  ®I/yMV  +  (alIvM^/dT,  (T  ”  To)  * 

CLVMV  "  Ww  +  1lYMVABEAB  +  lLYKMWV,  K  +  2LKMVWV,  K  *  (2,3) 

Equations  (2.1)  constitute  the  stress  equations  of  notion  and  charge  equation 
of  electrostatics  referred  to  the  position  coordinates  of  material  points 
before  the  thermally  induced  static  deformation  is  applied.  Which  are  called 
reference  coordinates  and  are  denoted  by  X^.  Equations  (2.2)  are  the  linear 
electroelastic  constitutive  equations  and  Eqs. (2.3)  contain  the  definitions 
associated  with  the  effective  elastic  coefficients  appearing  in  (2.2).  In 

A,  ~ 

(2.1)  -  (2.3)  K.  ,  A  and  u  denote  the  components  of  the  small  field  Piola- 
Uy  L  Y 

Kirchhoff  stress  tensor  which  is  asymmetric,  reference  electric  displacement 
vector,  and  mechanical  displacement  vector,  respectively;  p°  and  cp  denote 
the  reference  mass  density  and  small  field  electric  potential,  respectively, 
§LYMv  and  §LYMvAB  <3enote  t^ie  second  and  third  order  elastic  constants, 
respectively,  e_„.  and  e_„  denote  the  piezoelectric  and  dielectric  constants, 

JjMV  LM 

respectively,  and  Tq  and  T  denote  the  reference  and  present  temperature, 
respectively. 

In  this  description  the  present  position  £  of  material  points  is  re¬ 
lated  to  the  reference  position  X  by 

X(XL,t) -X  +  w(XL) +u(XL,t)  ,  (2.4) 

Where  w  denotes  the  displacement  due  to  the  thermally  induced  static  deforma- 
tion.  In  (2.3)  T*  and  E?_  denote  the  components  of  the  thermally  induced 

LiM  Ad 

static  biasing  stress  and  strain,  respectively,  and  for  small  strains  and 


changes  in  temperature  from  the  reference  temperature  T  ,  we  have 


4. 


BAB  “  2  (WA,B  +  WB,A)  » 


(2.5) 


Tm  “  §i*fsers  “  vlm(t  ~  To)  * 


where  vTU  denotes  the  thermoelastic  coupling  coefficients.  Which  are  related 

in 

to  the  usual  coefficients  of  linear  thermal  expansion  o'  by 

VI*=CW*JK-  (2*7) 

The  upper  cycle  notation  for  many  dynamic  variables  and  *-he  capital  Latin 
and  lower  case  Greek  index  notation  are  being  employed  for  consistency  with 
the  notation  of  Refs. 5  and  6,  as  is  the  remainder  of  the  notation  in  this 
section.  Hie  fact  that  the  capital  Latin  and  lower  case  Greek  indices  refer 
to  the  reference  and  intermediate  position  coordinates,  respectively,  is  not 
important  here  and  in  this  work  they  may  be  used  interchangeably.  We  employ 
Cartesian  tensor  notation,  the  summation  convention  for  repeated  tensor 
indices,  the  convention  that  a  comma  followed  by  a  Latin  index  denotes 
partial  differentiation  with  respect  to  reference  coordinates  and  the  dot 
notation  for  differentiation  with  respect  to  time.  Since  we  are  considering 
a  homogeneous  stress  free  biasing  state  resulting  from  a  homogeneous  temper¬ 
ature  change  (T  -  T  ),  from  (2.6)  and  (2.7)  we  obtain 
o 

E*s  '  «rs<t-V  •  ,2-8» 

4 

Since  it  has  been  shown  that  a  static  homogeneous  infinitesimal  rigid  rota¬ 
tion  has  no  influence  on  any  results,  without  any  loss  in  generality  we  may 
choose  the  homogeneous  rigid  rotation  to  vanish  and  substitute  from  (2.8) 


into  (2.3) j  to  obtain 


LyMV  “  (§LYMvABaAB  +  2LYKMQ'vK+  ^LKMvS'k'  (T  "  To) 


5 


3.  Perturbation  Equations 

For  purely  elastic  nonlinearities,  Which  are  the  only  ones  of  interest 

here,  the  equation  for  the  first  perturbation  of  the  eigenvalue  obtained  from 

0 

the  perturbation  analysis  mentioned  in  the  Introduction  may  be  written  in 


the  form 


AV/VX  *  1^/25  VL  ,  V  =  VL  +  AV , 


where  and  V  are  the  unperturbed  abd  perturbed  surface  wave  velocities  at 
T  and  T,  respectively,  ?  is  the  unperturbed  propagation  wavenumber  and 


Hl-+Ji?VSY,t.aV’ 


where  R  is  the  volume  enclosed  in  a  wavelength  in  a  unit  length  of  the  sur¬ 
face  wave  from  the  top  of  the  substrate  down.  In  (3.2)  denotes  the  nor¬ 
malized  mechanical  displacement  vector,  and  denotes  the  portion  of  the 
Piola-Kirchhoff  stress  tensor  resulting  from  the  biasing  state  and  the  change 
in  the  elastic  constants  with  temperature  in  the  presence  of  the  g^. 


and  is  given  by 


where 


*LY  =  ^LYMor +  ^2LYMc^  gar,  M  ’ 


V**=(d2^/aT)(T'V>  (3-4 

and  the  dc^j^dT  are  the  first  temperature  derivatives  of  the  fundamental 
elastic  constants  of  quartz,  which  are  tabulated  in  Ref. 4.  The  normalized 
eigensolution  g  and  f  is  defined  by 


*v"  >  n2  =  J  p\VV> 


6. 


f  • 


where  and  cp  are  the  mechanical  displacement  and  electric  potential, 

respectively,  which  satisfy  the  equations  of  linear  piezoelectricity,  i.e., 

Eqs.  (2.1)  and  (2.2),  with  (2.3)  and  (2.9),  at  T=T  subject  to  the  boundary 

o 

conditions 


■  Av  '  0  -  ‘ Vl  '  *„£  -  “  Vl  _  0  > 


(3.6) 


at  T  =  TQ,the  first  of  which,  with  (2.2)1,  (2.3)  and  (2.9),  takes  the  form 

V!lyMVUv,M+W,«>-°-  (3.7) 


The  second  equation  in  (3.6)  is  a  consequence  of  Laplace's  equation  in  free 
space,  the  two  electrical  continuity  conditions  at  the  surface  of  the  sub¬ 
strate  and  the  form  of  the  surface  wave  solution  appearing  in  (4.1).  In 

(3.6)  and  (3.7)  N  denotes  the  unit  normal  to  the  free  surface  of  the  sub- 
L 

strate  at  T=T  and  e  is  the  permittivity  of  free  space, 
o  o 


4.  Temperature  Dependence  of  Surface  Wave  Velocity 

A  schematic  diagram  of  the  free  surface  of  the  substrate  at  the  refer¬ 
ence  temperature  T  =  Tq  is  shown  in  Fig. 1.  The  solution  for  surface  waves 
propagating  in  the  X^-direction  may  be  written  in  the  form 


r  ,  ,  ,  .  ,  ,  iP  ?X  i?  X  -V.t) 

\  _  (m)  (m)  (m)  nr  v  *  t  1 

(UY’CP)  L  C  (\  ,B  )e  e  » 

m=l 


(4.1) 


where  the  unperturbed  surface  wave  velocity  v  at  T  =  Tq  is  determined  numer¬ 
ically  by  trial  and  error  so  that  (2.1)  and  (3.6),  with  (2.2),  (2.3)  and  (2.9), 
10-12 

T  =  Tq  are  satisfied  .  When  V  is  obtained  the  equations  determine  the 
Pm>  A^,  and  C^.  Since  the  straight-crested  surface  wave  solution 

depends  only  on  the  two  coordinates  and  Xv,  for  integration  over  a  wave¬ 


at 


length  Bq. (3.2)  takes  the  somewhat  reduced  form 


7. 


TT/| 

H  =  -  f  dX  [  dX  (V1  g  +Kng  +  K  g 
1  J  V  J  tL  TT  T,T  TV  V,T  Tff  U,T 

"  -tt/5 

•vn  '■v-  j 

k  ,  g  +  K  g  +Kg  , 
vt  t,v  vv  v,  v  vcr  a,  vj  ’ 


(4.2) 


where  X  denotes  the  axis  normal  to  X  and  X  and  in  (4.2)  we  have  introduced 
CT  TV 

the  convention  that  repeated  Greek  indices  are  not  to  be  summed.  For  this 
case  of  surface  waves  the  normalization  integral  in  (3.5)3  takes  the  form 


2  o 
N  =  p 


J 


dX 


tt/I 

I.  V 


-n/S 


V  VdXT  = 


p°TTi  «  «  c(m>^m)C(n)  ^ 


l  l} t 

m=l  n=l 


<0  -P*) 

m  n 


(4.3) 


where  the  *  denotes  complex  conjugate.  Hie  substitution  of  (4.1)  and  (4.3)  into 

(3. 5)^  which,  with  (2.9)  and  (3.4),  is  then  substituted  into  (3.3),  which  is 

then  substituted  along  with  (3. 5) ^  into  (4.2),  the  integration  of  which 

yields  the  expression  for  H^,  which  enables  the  change  in  surface  wave 

velocity  to  be  calculated  from  (3.1).  The  resulting  expression  for  is 

extremely  lengthy  and  not  terribly  revealing.  Consequently,  we  do  not 

bother  to  present  it  here.  Nevertheless,  it  is  clear  from  (4.2)  and  (3.3), 

with  (3.4)  and  (2.9),  along  with  (3.5)^,  (4.1)  and  (4,3)  that  the  resulting 

13  14 

express.' on  for  is  linear  in  (T  -  Tq)  ,  depends  on  the  second  and  third 

order  elastic  constants,  the  coefficients  of  linear  expansion^5  and  the  recently 
4 

determined  first  temperature  derivatives  of  the  fundamental  elastic  constants 

of  quartz,  as  well  as  on  the  known  surface  wave  solution  at  T=T  which 

o7 

13 

depends  on  the  linear  piezoelectric  constants  . 

The  change  in  velocity  AV  determined  from  (3.1)  is  the  change  in 
16 

natural  velocity,  which  is  related  to  the  measured  unchanged  geometry  at 
T*Tq.  Since  it  is  not  feasible  to  continually  measure  the  changing  lengths 


with  temperature  and  time  intervals  are  usually  measured  rather  than  the 


8 


actual  velocity  of  the  wave,  the  natural  velocity^  is  the  most  purposeful 
one  to  define.  However,  since  previous  workers  have  employed  the  change  in 
actual  velocity  calculated  from  the  measurements,  we  will  calculate  the 
change  in  actual  velocity  from  the  change  in  natural  velocity.  The  relation 
between  them  is^ 

AV/V.  •  AV/V.  +  eV.  ,  t  no  sum,  (4.4) 

a  1  ±  tt 

where  the  repeated  Greek  indices  in  (4.4)  are  not  to  be  sunned  and  AV  is 

A 

the  change  in  actual  velocity  as  defined  in  the  Appendix.  Since  reference 
coordinates  and  natural  velocity  are  employed  here,  the  geometry  remains 
fixed  and  the  temperature  coefficient  of  delay  (TCD)  is  given  by 

TCD*At/t(T-T  )  -  -  AV/V,  (T  -  T  ),  (4.5) 

o  i  o 

where  t  is  the  delay  time.  Equation  (4.5)  shows  that  when  reference  coordi¬ 
nates  and  natural  velocity  are  employed,  the  temperature  coefficient  of  delay 
is  simply  the  negative  of  the  temperature  coefficient  of  velocity.  Substi¬ 
tuting  from  (4.4)  into  (4.5)  and  employing  (2.8),  we  obtain 

TCD*  -  AV  /V.  (T  -  T  )  +a  ,  t  no  sum,  (4.6) 

a  l  o  tt  7 

which  is  the  well-known  usual  relation  between  the  TCD  and  the  temperature 

coefficient  of  actual  velocity  referred  to  the  variable  temperature  dependent 

intermediate  configuration  of  the  crystal. 

Calculations  have  been  performed  for  various  propagation  directions  on 

various  orientations  of  quartz  substrates  and  the  results  are  plotted  in 

Figs.2-8.  Figures  2  and  3  show  the  TCD  and  AV  /V.  (T  -  T  )  as  a  function  of 

a  i  o 

propagation  direction  for  surface  waves  on  AT-cut  and  AC-cut  quartz  substrates, 
respectively,  at  25° C.  Ihe  dotted  curve  shows  the  average  of  the  calculated 
values  at  0°C  and  50°C  from  Ref. 1  and  the  circles  are  the  average  of  the 
experimental  values  at  0°C  and  50° C  also  from  Ref.l.  It  is  clear  from 


Figs. 2  and  3  that  the  calculations  performed  here  using  the  proper  nonlinear ly 
based  formalism  are  in  substantially  better  agreement  with  the  experimental 
data  of  Ref. 1  than  the  calculations  performed  in  Ref.l  using  an  incomplete 
linearly  based  description.  Figure  4  shows  the  calculated  TCD  for  surface 
waves  propagating  along  the  digonal  axis  of  rotated  Y-cut  quartz  substrates 
as  a  function  of  rotation  angle  at  25°c.  the  dotted  curve  is  obtained  from 
the  calculations  of  Ref. 1  in  the  same  manner  as  in  Figs. 2  and  3.  Figures  5 
and  6  show  the  calculated  TCD  as  a  function  of  propagation  direction  for 
surface  waves  on  X-cut  and  Y-cut  quartz  substrates,  respectively,  at  25°  C. 

The  dotted  curves  are  obtained  from  the  calculations  of  Ref.l  in  the  same 
manner  as  in  Figs. 2  and  3.  Figures  7  and  8  show  the  temperature  coefficient 
of  delay  and  associated  power  flow  angle  as  a  function  of  propagation  direc¬ 
tion  for  two  particularly  interesting  orientations  of  doubly-rotated  quartz 
substrates,  namely  those  having18  tp*30°,  0  »  34.2^  and  cp-5°,  0  «  47.7° , 
respectively. 
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FIGURE  CAPTIONS 


Schematic  Diagram  Showing  the  Free  Surface  of  a  Semi- 
Infinite  Solid 

Temperature  Coefficients  of  Actual  Velocity  and  Delay  for 
Surface  Waves  on  AT-Cut  Quartz  as  a  Function  of  Propagation 
Direction  Relative  to  the  Digonal  Axis  at  25°  C.  The  dotted 
curve  shows  the  average  of  the  calculated  values  at  0°Cand 
50°  C  from  Ref.l.  The  circles  are  the  average  of  the  experi¬ 
mental  values  at  0°C  and  50° C  from  Ref.l 

Temperature  Coefficient  of  Actual  Velocity  and  Delay  for 
Surface  waves  on  AC-Cut  Quartz  as  a  Function  of  Propagation 
Direction  Relative  to  the  Digonal  Axis  at  25° C.  The  notation 
convention  is  the  same  as  in  Fig. 2. 

Temperature  Coefficient  of  Delay  for  Surface  Waves  Propa¬ 
gating  along  the  Digonal  Axis  of  Rotated  Y-Cuts  of  Quartz 
as  a  Function  of  Rotation  Angle  at  25°C.  The  notation  convention 
is  the  same  as  in  Fig. 2. 

Temperature  Coefficient  of  Delay  for  Surface  tteves  on  X-Cut 
Quartz  as  a  Function  of  Propagation  Direction  Relative  to 
the  Y-Axis  at  25° C.  The  notation  convention  is  the  same 
as  in  Fig. 2. 

Temperature  Coefficient  of  Delay  for  Surface  waves  on  Y-Cut 
Quartz  as  a  Function  of  Propagation  Direction  Relative  to 
the  X-Axis  at  25° C.  The  notation  convention  is  the  same 
as  in  Fig. 2. 

Tempera ture  Coefficient  of  Delay  and  Power  Flow  Angle  as  a 
Function  of  Propagation  Direction  Relative  to  the  Axis  of 
Second  Rotation  at  25° C  for  the  Doubly-Rotated  Cut  of 
Quartz  Having  cp=5°,  8=47.7°. 

Temperature  Coefficient  of  Delay  and  Power  Flow  Angle  as 
a  Function  of  Propagation  Direction  Relative  to  the  Axis 
of  Second  Rotation  at  25° C  for  the  Doubly-Rotated  Cut  of 
Quartz  Having  ep=  30°,  9  =  34.25°  . 

Schematic  Diagram  Showing  the  Relation  Between  the  Surface 
Normal  and  Surface  wave  Propagation  Direction  and  wave¬ 
length  in  the  Undeformed  Reference  State  at  Temperature  T 
and  the  Deformed  State  at  Temperature  T. 
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APPENDIX 


In  this  Appendix  we  derive  the  relation  between  the  change  in  actual 

and  change  in  natural  velocity  for  surface  waves,  the  procedure  employed 

19 

follows  that  of  Thurston  who  treated  plane  waves,  the  treatment  of  the 

surface  wave  case  is  a  bit  more  cumbersome  because  of  the  existence  of  the 

free  surface.  Consequently,  in  this  Appendix  the  distinction  between  capital 

Latin  and  lower  case  Greek  indices  is  important. 

Let  V.,  and  cr,  be  unit  vectors  normal  to  the  surface  and  in  the  direction 
V  V 

of  the  phase  fronts  of  the  surface  wave,  respectively,  both  at  temperature  T, 
as  shown  in  Fig. 9.  Clearly,  then  the  propagation  direction  is  given  by 


TT  *>  VXo 

'V  'V 


f 


TT 


SopVp ' 


(Al) 


where  e^p  is  a  skew-symmetric  tensor.  In  the  undeformed  state  at  tempera¬ 
ture  Tq  the  unit  vectors  v  and  a  were  in  the  directions  N  and  S,  respectively, 
in  Fig. 9.  We  denote  the  vector  orthogonal  to  N  and  S  by  P,  which  is  in  the 
natural  propagation  direction,  and  we  write 


P-NXS,  P=E  NS 
~  ~  ~  ’  X  KIM  L  M 


(A2) 


where  E  is  a  skew-symmetric  tensor.  In  the  deformation  accompanying  the 

lvI/4 

temperature  change  from  Tq  to  T  the  line  element  in  the  direction  P  at  Tq 

takes  the  direction  ji  at  T,  which  in  general  differs  from  the  direction  tt. 

The  intermediate  directions  v,  a  and  u  may  be  expressed  in  terms  of  the 

20 

reference  directions  N,  S  and  P  by  means  of  the  relations 

~  cw 


V  CT  „  ^Y.I.PL 

y  B(io  *  °y  c(s)  '  N  c(p)  ’ 


(A3) 


where 


A2 


B  "  "AcCKMNM  *  C  "  ^SKCKMSM  *  C  “  ,,^>KCKXPM  *  (A4) 


and 


ckm“  ^a,K^a,M»  CKM“XK,orXM,o  » 


(A5) 


where  “  ?y (XR)  is  the  intermediate  coordinate  at  T  and  is  related  to  the 
quantities  appearing  in  (2.4)  by 


V  WW  - 


(M) 


where  6^L  is  a  Kronecker  delta,  which  is  required  for  consistency  with  the 
notation  employed5,  and  a  comma  followed  by  a  Greek  index,  say  P,  denotes 
partial  differentiation  with  respect  to  the  intermediate  coordinate 

P 

Substituting  from  (A3)  into  (Al),  we  obtain 


(A7) 


which  is  clearly  different  from  (A3)3> 

Since  the  natural  (or  reference)  wavelength  L  is  in  the  direction  P, 
20 


from  deformation  theory  and  (A3) we  have 

where  i  is  the  deformed  length  that  was  L  when  undeformed  and 

l  -  C(P)L. 


(A8) 


(A9) 


Inasmuch  as  L  is  the  undeformed  (or  reference)  distance  between  two  phase 
fronts,  l  is  the  deformed  line  between  the  same  two  (deformed)  phase  fronts, 
which,  however,  is  not  normal  to  the  phase  fronts  because 

u  *  a?  0  .  (A10) 

K*  iv 

Since  the  distance  between  two  deformed  phase  fronts  normal  to  the  phase 
fronts  is  the  actual  wavelength  X,  the  component  of  ij£  in  the  direction  tf  is 


*  P 


A3. 
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the  actual  wavelength  ,  and  we  have 

X  -  Jtjj  •  tt  .  (All) 

Substituting  from  (A7)  and  (A8)  into  (All)  and  rearranging  terms,  we  obtain 

(A12) 

.  ij  jj,  or 

20 


Utilizing  the  well-known  relations 


'60Y 


^6,R^,K^V,M*  J  ERKM»  ^6,RXR,a”  6cr6  , 


we  obtain 


sctPY53,K^Y,m”  J\ oERKM  ’ 


where,  of  course 

a  det  I  . 

Substituting  from  (A14)  into  (A12)  and  rearranging  terms,  we  obtain 

k“ulWKVR,A,aVB'",c'S)  ■ 

which,  with 

SXP-N.  ESP  « N  , 

~  ~  ~  *  RKM  KM  R 7 

and  (A4)  and  <A5)2,  yields 

x-lJVn)/c(S)> 


(A13) 


(A14) 


(A15) 


(A16) 


(A17) 


(A18) 


which  holds  for  arbitrarily  large  deformations.  Before  reducing  to  infinites¬ 
imal  deformations  we  record  the  well-known  relation 


EKM  ”  2  (CKM  ”  5KM)  * 


(A19) 


When  the  deformation  is  infinitesimal  and  terms  higher  than  linear  in 
w  may  be  ignored,  from  (A5),  <A6)  and  (A19)  we  obtain 

U,R 


A4 


5  <wit,H*’V)  • 


(A20) 


'KM 


6km+2ekm> 


. 

KM 


6  KM  ”  2EKM  * 


(A21) 


Substituting  from  (A21)  into  (A4)‘,  utilizing  the  fact  that  N,  S  and  P  are 
unit  vectors  and  making  a  Taylor  expansion  of  the  radicals  while  retaining 
only  linear  terms,  we  obtain 


B 


(N) 


1  -  E 


(N)  1 


,<S) 


1+E 


(S)l 


.CP) 


1  +  E 


(P) 


(A22) 


where 


,  (N)  1 


VkmV 


,(S)1 


VtoAi' 


E 


(P)l 


pkekmpm ; 


(A23) 


and  E(N)1,  E  ^ 1  and  E^1  are  the  extensional  strains  in  the  directions 
N,  S  and  P,  respectively.  Similarly,  substituting  from  (A6)  into  (A15), 
expanding,  retaining  only  linear  terms  in  w  and  employing  (A20),  the 

K 

invariance  of  the  trace  of  the  strain  tensor  and  (A23) ,  we  obtain 


J 


1 


i  +  e(n)1  +  e(S)1  +  e(p)1. 


(A24) 


Finally,  substituting  from  (A22)  and  (A24)  into  (A18),  making  a  Taylor  ex- 
(S)  1 

pansion  in  the  E  ,  multiplying  out  and  retaining  only  linear  terms  in 
(• )  1 

the  E  ,  we  obtain 

(P)l 

X-  L(l  +  E'  '  )  ,  (A25) 


which  is  the  equation  that  has  always  been  employed  without  a  complete 

demonstration  of  its  validity.  Since  the  frequency  is  the  same  whether  the 

reference  coordinates  X  at  T  or  the  deformed  coordinates  E  at  T  are 

~  o 

employed,  we  have 

V  A  -  V/L  , 

s 


(A26) 


A5 


where  V  is  the  actual  velocity  and  V  is  the  natural  velocity  of  the  surface 

A 

wave.  Since 

V  -  V  +fcV  ,  V-V.+AV,  (A27) 

a  X  A  X 


where  V  is  the  unperturbed  velocity  at  T  and  both  AV  and  AV  are  small 
J.  o  * 

compared  with  V^,  from  (A25)  -  (A27),  we  have 


AV 


AV  (P)l 

v“  > 

1 


(A28) 


(P)  1 

because  E  «  1.  Equation  (A28)  is  Eq. (4.4). 

Since  we  have  presented  the  treatment  in  this  Appendix  leading  to 
Eq. (A18) ,  from  which  we  have  obtained  the  reduced  form  (A25)  for  the  case 
of  infinitesimal  strain,  for  completeness  we  record  without  proof  the 
limiting  forms  of  v,  a,  and  tt  for  the  case  of  infinitesimal  strain.  When 
we  ignore  the  distinction  between  lower  case  Greek  and  capital  Latin  indices 
these  forms  are; 


ttl”  pl“  slsmwk,  mpk  +  Wm,  k^k  * 


(A29) 


